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Abstract 



We have investigated the classical stability of magnetically charged black 
p-brane solutions for string theories that include the case studied by Gregory 
and Laflamme. It turns out that the stability behaves very differently depend- 
ing on a coupling parameter between dilaton and gauge fields. In the case of 
Gregory and Laflamme, it has been known that the black brane instability 
decreases monotonically as the charge of black branes increases and finally 
disappears at the extremal point. For more general cases we found that, 
when the coupling parameter is small, black brane solutions become stable 
even before reaching to the extremal point. On the other hand, when the cou- 
pling parameter is large, black branes are always unstable and moreover the 
instability does not continue to decrease, but starts to increase again as they 
approach to the extremal point. However all extremal black branes are shown 
to be stable even in this case. It has also been shown that main features of 
the classical stability are in good agreement with the local thermodynamic 
behavior of the corresponding black hole system through the Gubser-Mitra 
conjecture. Some implications of our results are also discussed. 
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I. INTRODUCTION 



Black holes are very important objects for studying and understanding both general 
relativity and string theory. Brane configurations in string theory have been studied from 
various aspects, and play important roles in the context of gauge/gravity duality From 
the AdS/CFT correspondence it has been also discussed that the IR limit of non supersym- 
metric Yang Mills theories corresponds to black hole configurations in AdS space which can 
be regarded as the near horizon limit of non extremal D-brane configurations . From the 
study of Yang Mills theories we can reveal properties of black holes, and it is discussed that 
the absence of tachyonic glueball implies the stability of corresponding black hole config- 
urations 0. On the other hand, recently the issue about the black string instability and 
its evolution has been of great interest. It was widely believed that the evolution of the 
instability in a black string results in bifurcations of the horizon finally. However, Horowitz 
and Maeda have shown that the horizon of a black string cannot pinch off in a finite affine 
time [Q. Therefore an inhomogeneous black string may be the final stable configuration. 
Recently some authors discuss about the properties of inhomogeneous black strings and their 
existence 

The Schwarzschild black hole is known to be stable under linearized perturbations in 
general relativity [^. However, Gregory and Lafiamme discussed the stability of black 
branes and found that a foliation of Schwarzschild black holes (Sch. x S") is unstable if 
the compactification scale (S"" radius) is larger than the order of the Schwarzschild radius 
- the so-called Gregory-Lafiamme instability One may, roughly speaking, expect that 
charge prevents a black brane from being unstable because of some repulsive forces among 
charges as the event horizon shrinks. However they also showed that some charged black 
branes in ten dimensional spacetime always have instability modes all the way down to 
the extremal limit []T0[. So, it is unlikely that the presence of charge is possible to remove 
the instability in black branes they considered. On the other hand, extremal black branes 



are shown to be stable |TT|. The stability of black branes in anti-de Sitter spacetime has 
also been investigated. In five-dimensional anti-de Sitter spacetime, black strings that are 
foliations of four-dimensional Schwarzschild black holes IT^ or de Sitter-Schwarzschild black 



holes [0 are recently shown to be unstable as well. However the black string that is a 
foliation of four dimensional anti-de Sitter-Schwarzschild black holes becomes stable as the 



horizon radius is larger than the order of the AdS4 radius ||T^. BTZ black strings in four 



dimensional spacetime turn out to be stable always |jT4|. 

In string theory context, black branes that Gregory and Lafiamme considered p!0|JTl[] are 
those having magnetic charge with respect to Neveu-Schwarz gauge fields {i.e., the case of 
a = —2 in Eq. (0) below). Thus, for example, black p-branes carrying charge with respect 
to Ramond-Ramond gauge fields {i.e., the case of a = 0) are not covered in their work. 
Properties of D-branes and their non extremal extensions have many applications and are 
important to study. Charged black brane solutions Gregory and Lafiamme considered might 
not be general enough to show other possible interesting stability behavior. Therefore, it 
will be of interest to see whether or not the stability behavior drastically changes when 
more general types of charged black brane solutions are considered. In fact we find different 
behaviors depending on the coupling between dilaton and gauge fields. For instance charged 
black branes become stable even far from the extremal point provided that this coupling 
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is smaller than a certain critical value. In this paper, we study the classical stability of 
a wider class of charged black brane backgrounds which include various black p-branes in 
string theory as special cases. 

Another interesting feature of black hole systems is that they can be regarded as ther- 
modynamic systems. A black hole has entropy, temperature, and other thermodynamic 
quantities. The area theorem states that the total area of black holes cannot decrease, i.e., 
6 A > 0, and it might be interpreted as the total entropy cannot decrease 55* > in the black 
hole thermodynamics. The behavior of instability explained above has been commonly con- 
sidered in the viewpoint of thermodynamics as a reflection that a black string is entropically 
less-favorable compared with multi black holes having the same mass and conserved charges. 
However, the classical area theorem does not necessarily require the transition of a black 
brane into a configuration having larger horizon area. Thus, this global entropy argument 
is too naive and cannot be applied to all cases. For example, one can show that a black 
string in five-dimensional AdS space can be stable even if its entropy is smaller than that of 



five-dimensional AdS black holes [13|. As pointed out by Reall [15|, the relevant one should 



be the local thermodynamic behavior. 

Recently, Gubser and Mitra gave such refinement of the entropy argument and conjec- 
tured that a black brane with a non-compact translational symmetry is classically stable 



if and only if it is locally thermodynamically stable (GM conjecture) |16]. Reall argued 
general correspondence between the classical instability and the presence of a negative eigen 
mode in the Euclidean Einstein-Hilbert action [|15[. The negative eigen mode gives an 



imaginary contribution to the path integral and then it is closely related to the local ther- 
modynamic instability. In fact the local thermodynamic instability implies the existence of 
the negative mode(s). However, there is no rigorous proof for the converse: the existence 



of a negative mode indicates the thermodynamic instability. For Schwarzschild |TB| and 



Reissner-Nordstrom |jT6| black holes in AdS space and a Schwarzschild black brane enclosed 
in a finite cavity |T^, the equivalence between non-existence of negative modes and local 
thermodynamic stability are already examined by performing both classical and thermody- 
namic stability analyses explicitly. 

However, it is worthwhile to study classical stability of black branes in the connection 
with black hole thermodynamics in a more general context. Charged black brane solutions in 
type II supergravity are particularly interesting. As pointed out by Reall magnetically 



charged p-branes with p < A (DO, Fl, Dl, D2, D4) have sign changes in the specific heat 
as the charge increases. Although in Ref. []T3[ the essential feature of local thermodynamic 
stability for those black branes has been analyzed and a sketch of the proof for the GM 
conjecture is given in the context Reall considered, the actual analysis for the classical 
stability has not been performed. So, it has not been examined explicitly as yet how well 
classical and thermodynamic stabilities actually agree with each other. This is another main 
motivation of our study. By performing the classical stability analysis, we give an explicit 
check for the validity of the GM conjecture. As shall be shown below, the agreement between 
our numerical results for classical perturbations and the thermodynamic stability through 
the GM conjecture is very good. We also show some features in the classical stability 
behavior that are not easily predicted by the thermodynamic behavior only. 

This paper is organized as follows. In section |I|, we summarize the thermodynamic 
stability behavior of black branes we are considering. Section fT| is devoted to the numerical 
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analysis for the classical stability. We give conclusions and discussions in the last section. 



II. LOCAL THERMODYNAMIC STABILITY OF BLACK P-BRANES 



In general, it is not easy to check whether or not black string or brane background 
spacetimes are stable under classical linearized perturbations. Such classical analysis in- 
volves numerical calculations for most cases. As stated above briefly, however, there has 
been suggested recently a simple method from which one can predict about the basic be- 



havior of the classical stability. This is the so-called Gubser-Mitra (GM) conjecture [jlG 



In this section, before performing the classical analysis in detail, we first consider the local 
thermodynamic behavior of a finite segment of a black p-brane background spacetime. 
We consider black p-brane solutions of the following action 



(1) 



Here F„ denotes a n-form field strength. /3, 7, and a are assumed to be arbitrary constants 
with 7 = 7 + (D — l)/3^/{D — 2) > 0. However, by taking the conformal transformationB 
gMN = gMN and then rescaling the dilaton field (p = 0/^/27, this action can be 

written as 



I 



R 



2n\ 



where 



a+ {D -2n)p/{D -2) 



(2) 



(3) 



Note that the sign of the parameter a is not relevant here since the action is invariant 
under a — >■ —a and — > —(p. In the following analysis, we will consider only magnetically 
charged black brane solutions of this action, not electrically charged ones. It is because 
electrically charged solutions can be obtained from magnetically charged ones by dualizing 
the n-form field F. Moreover, the magnetically charged case is technically easier to treat in 
the perturbation analysis.i The metric of magnetically charged black p-brane solutions is 
given by Em 



Note that the classical stability of black branes does not change under a conformal transfor- 
mation as long as the conformal factor is regular. Moreover, black hole thermodynamic quanti- 
ties such as temperature |^^, entropy, and specific heat are invariant under stationary conformal 
transformations . 

^Gauge perturbations can be set to be zero for the sphere directions, and are decoupled from 
metric and scalar perturbations for other directions. 
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ds"^ = — f 1 H = sinh^ u 



( k 
+ 1 H — = sinn u 

V r'^ J 
e 27r<P = 1 + ^ sinh /i, 



"A(D-2) ( k ^ 

Udt^ + I 1 + — sinh^ 



A(D-2) /^^2 

[it 



Ad 



A(D-2) 



6ijdz''dz\ 



U 



k 



A = a^ + 



(4) 



The n-form field strength is proportional to the volume form on S^. Here d = n — 1, 
d = p+l, D = d + d + 2 = 2 + n + p, and the coordinates are {x^} = {x'^, z^} = {t, r, x™, z^} 
with m = 1, ■ ■ ■ ,n and i = 1, ■ ■ ■ ,p. Thus, the coordinates {z*} describe the p-dimensional 
spatial worldvolume, and {x™} the n-sphere. The spatial worldvolume is not compactified 
in our consideration. 

Note that the metric components for the p-dimensional spatial worldvolume directions 
above have an overall multiplication factor depending on the r-coordinate. This multipli- 
cation factor makes the spatial worldvolume non-flat, and turns out to give some compli- 
cations in the perturbation analysis. As pointed out by Reall |15], such complications can 
be easily avoided by performing an appropriate conformal transformation again such as 
Qmn = e^^"~^^'^^^^~'^^"'gMN- Then the action becomes 



where the constants are now functions of a single parameter a as follows: 

n-1 1 ^ {D-2n){n-l] 

2 ~ 
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a = a 



a 2 {D-2)a^ ' {D - 2)a 

And the metric for the black p-brane background can be written as 

ds"^ = -Udt^ + V-^dr^ + R^dnl + 6i4z'dz\ 



where 



_i_ (l + ;^sinhV^ 



k/r" 



= h + ^sinhV 



A 



(5) 



(6) 



(7) 



(8) 



The mass M and magnetic charge Q per unit p-volume are given by ^2 



Ad 

M = k{d + l + — sinh^ fx) and Q 



dk 



sinh 2/i, 



(9) 



respectively. 

Note that the case of /3 = 2, 7 = —4, a = —2, and = 10 in Eq. corresponds to the 
action for which the classical stability of magnetically charged black p-brane solutions has 
been studied by Gregory and Lafiamme (GL) [|TU|,|TT|. Thus, as pointed out by Reall |]T3 



NS5-brane {i.e., = 3) is the only black brane of the type II supergravity covered by their 
work. This case corresponds to a = (1 — 'n,)/2 in the conformally transformed action in 
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Eq. d^). F-string by dualizing the Neveu-Schwarz gauge field can also be described by the 
action in Eq. with a different coupling to the dilaton, i.e., a = 2 and so a = (9 — n)/2. 
Hence F-string {i.e., n = 7) can be covered. When /3 = 2, 7 = —4, a = 0, and D = 10, it 
gives a = (5 — n)/2. Solutions in Eq. (0) for this case cover Dp-branes with p = 0, 1, 2, 4, 5, 6 
of the type II supergravity carrying magnetic charge. Therefore, our study contains a broad 
class of black brane solutions, including most of black branes in the type II supergravity. 

Now let us consider a black p-brane with unit worldvolume. Being regarded as a thermal 
system, it has entropy and temperature given by 

S ~ (cosh/i)^ and T = (cosh/i)"^ , (10) 

respectively. Here r„ = k^/'^ is the horizon radius. The extremal limit of the magnetically 
charged black brane solution is — > and fi 00 with keeping mass and charge in Eq. 
finite [i.e., ke^^ fixed). Note that in the extremal limit both temperature and entropy 
go to zero except for the case of A = 2d, i and the solution becomes the B.P.S. state 
which is known to be stable even at the quantum level. In order for this system to be 
stable thermodynamically, the entropy functional given above should be a local maximum 
in thermodynamic phase space. Defining the Hessian of the system as 

TT i dAP dMdQ \ /-,-,\ 
\ d&WQ / 

the maximum entropy is guaranteed provided that 



< and det (H) > 0, (12) 



Q 



which become equivalent to 



Here $^ denotes the magnetic potential energy at the horizon. Using Eqs. (P,[TO|), we obtain 
explicitly 



^ -4.4« (cosh ,)* ^"^j^;^'^:^"-"^'' , (14) 

2(3 + (A — 2dj cosh 2/i 
fd^A r-'^A cosh 2/i 

\dQ )^ 2rf+ (A-2(i)cosh2/i' 

Since the term A + d{A — 2) is positive definite, one can find that Cq and ((9$H/c^<5)r 
always have opposite signs in the above. Consequently, two conditions in Eq. (|TB|) cannot 



^For the case of A = 2d, interestingly the temperature does not vanish although the entropy still 
goes to zero in the extremal limit. 
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be satisfied simultaneously. In other words, the black brane system in the consideration 
cannot be thermodynamically stable for processes in which both mass M and charge Q 
vary. Therefore, according to the GM conjecture, this may indicate that the black brane 



background is always unstable under classical perturbations. As mentioned by Reall [[15 
however, the charge of the black brane cannot fluctuate since there is no matter field carrying 
charge in our theory. Thus, only thermodynamic processes through which the charge does 
not change are relevant. Accordingly, the local thermodynamic stability of a black brane is 
determined solely by the sign of the specific heat Cq above. 



The specific heat in Eq. (IT^ is always negative if A — 2(i > 0, equivalently, \a\ > Uc 



where the critical value of a depends on the form n of the gauge field, or the dimension p of 
the spatial worldvolume as 

n-l D-3-P ^^g^ 



For \a\ < Ocr, however, the specific heat is negative if < fi < ficr, but positive if /Zcr < < 
oo. At /i = ficT, it diverges, i.e., Cq =Foo as /i — ficr T 0. Here the critical value of the 
parameter is 

2 [2{D - 3 - p)2 - (D - 2)a2] ^ ^ 

Based on the GM conjecture, therefore, it can be expected under classical linearized pertur- 
bations that a black p-brane background is always unstable provided \a\ > a^r, independent 
of its charge. For \a\ < a^, however, it becomes stable for fi > fi^r- For instance, black 



p-branes considered by Gregory and Lafiamme ||IO| , pA| have a = (1 — n)/2 as explained 
above. This value is exactly the same as the critical value for D = 10 in Eq. (|16D, i.e., 
= a^j,, implying negative specific heat always. Thus it is expected that all black branes 
considered by them, including NS5-brane, are unstable, and this behavior is shown by them. 
For F-string, we have a = 1 < = 3. For Dp-brane solutions of the type II string theory, 
we have a = {p — 3)/2 and = (7 — p)/2. Therefore, one can easily see that the specific 
heat changes sign for magnetically charged DO, Fl, Dl, D2, D4 branes, but not for D5, 
D6 branes as mentioned by Reall [1^. The case of D3-branes will not be considered in this 
paper since the action for the self dual gauge field is not known. Moreover, even if we simply 
use the action in Eq. (|I]) and impose the self duality condition, we encounter some problems 
for analyzing linearized perturbation equations. Finally we would like to point out that, 
except for cases of \a\ = Ocr, the specific heat goes to zero in the extremal limit, satisfying 
a formulation of the ordinary thermodynamic third law. 



III. CLASSICAL PERTURBATION ANALYSIS 



So far, we have analyzed the local thermodynamic stability of black p-brane systems that 
gives some hints at the stability of such background spacetimes through the GM conjecture. 
Now let us perform the classical stability analysis explicitly. We consider small metric 
perturbations about black p-brane background spacetimes and see whether or not there 
exists any mode that is regular spatially, but grows exponentially in time. 
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Since the background black p-brane spacetime in Eq. is independent of the time 
coordinate t and spatial worldvolume coordinates z^, one can assume that 



6(P = e^i+im.z^ f{r^ x"*), 5F = e™+*™'"'5F(r, x"') (18) 
for unstable mode solutions. Here > and the Kaluza-Klein (KK) mass rrii is a continuous 



real number. As explained in Refs. |ll^ . [15| , one can further use the diffeomorphism symmetry 



so that the scalar and vector parts of metric fluctuations Hmn are set to be zero for massive 



modes in the viewpoint of the Kaluza-Klein reduction |24 



H^i = Hi, = 0. (19) 

Since in general non-s-wave fluctuations are expected to be more stable than s-wave ones, 
we consider only spherically symmetric perturbations.^ Thus, as shown in Refs. the 
following gauge choice can be made 

Htm = Hrm = 0, = (20) 

for such perturbations. Other components Hu, Htr, H^r and the dilaton and field strength 
perturbations are also independent of the angular coordinates x^. That is, they are functions 
of the radial coordinate r only. 

For these forms of perturbations above, one can easily show that linearized equations for 
gauge fluctuations 5F are decoupled, i.e., VJv(e"'^5F^^i--P"-i) = g. Moreover, as in Ref . HTUI 



it can be shown that for background spacetimes in Eq. there exists no unstable mode 
solution for the spherically symmetric fluctuations. One can first show that the fluctuations 
proportional to the volume form on S"" should vanish if f2 ^ 0. Thus only remaining 
components 5Ftii...j„_i, 5-Frii-- j„_i, SFtrii--i„-2: and of such fluctuations could be non- 

zero. Here ii, ■ ■ ■ ,in-i are the spatial worldvolume coordinates z. By using the Bianchi 
identity V[Af'^-^Pi --Pn] = ^^r {NPi ■ ■ ■ Pn) = (trii ■ ■ -in-i) and the perturbation equations 
for (Pi ■ ■ ■ Pn-i) = {rii ■ ■ ■in-2), {ih ■ ■ ■'in-2) above, one obtain coupled equations for = 
J2imi6Ftii^...i^_2, fr = J2imi6Frii^...i^_2, and 5Ftri^...i^_^ that can be easily decoupled as 



i^f;' + (---)f;-U2 + - f, = o, (21) 




where m? = J^i"^!- Since both U and V are positive outside the horizon, one can see that 
the only regular solution of this equation is ft = 0, giving f,. = 5Ftri^...i^_^ = correspond- 
ingly. Similarly, one also finds that there are no non-vanishing regular solutions for other 
components of the fluctuations. Therefore we set 6F = in the following analysis. 
As shown in Ref. [|T5|, now the dilaton perturbation equation is given by 



^We will give an argument below that the consideration of s-wave perturbations only might be 
enough to show stability. 
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2 n - 11 



n 



(22) 



And linearized equations for metric perturbations H^i, are given by 

\j^H^, - 2 V(M V'^.)p + Vp V. - '^RpipK) + '^Rf.pu.HP'' 

+(3 (2 V(M - V'^mO Vp - 2/3 Vp V./ + 4(7 + f3^) V(p V.) / 



(n-1)! 



[n - 1)/7^'^F^,,,...,„_,F,.'^^ 



Al - An 



Ai ■••A„,_i 



pv 



The {yU^} and {ij} components are 



2(7 + /5')/V^0 

^; - 2/5/ 



0, 
0, 



(23) 



(24) 
(25) 



respectively. Here the covariant derivative and curvature tensors in the above equations are 
constructed by the metric g^j^y only. In other words, perturbation equations on D-dimensional 
black p-brane backgrounds become coupled second order differential equations for dilaton 
and gravitational fields with mass m on (D — ]9)-dimensional black hole backgrounds, which 
is expected from the viewpoint of Kaluza-Klein dimensional reductions. 

As pointed out in Ref. already, notice also that the dilaton perturbation equation 
Eq. (0) is actually not independent, but can be obtained from the trace of Eq. (^) and 
Eqs. ( |24| , p5|) . In Eq. (pSf) , one finds that the dilaton fluctuation / can be computed once the 



metric fluctuations H^j,y are obtained. Thus, we need to consider only the metric perturbation 
variables, which are Hu-, Htr, Hrr, and K. For the stability analysis, however, it suffices 
to consider the behavior of the so-called threshold modes, i.e., KK massive modes having 
= 0. This is because infinitely heavy massive modes cannot be expected to give instability. 
There should be a maximum value m* for the instability to exist. For unstable modes in the 
case of black brane systems, if they exist, it is known that the corresponding "imaginary" 
frequency fl usually starts to increase as the KK mass increases, reaches to some maximum 
value, but finally decreases to zero at the so-called threshold mass m*. Thus, as long as the 
set of equations in Eqs. ( p^ - |25|) allows a threshold mode solution with non-vanishing mass, 
the instability exists. When = 0, there is an additional merit that a further gauge choice 
is possible to set Htr = as in Ref. |]15|. Therefore, one finally has 



{Hl^) = diag (<^(r), V'(r), x(r), ■ ■ ■ , x(r)) • (26) 
That is, the perturbed metric for threshold modes can be expressed in the following form 
ds^ = -U(l + ipe'"^''') dt" + V-^ (l + ^e''"'"") dr^ + R^l + xe'""^"") d^l + rfzl (27) 
Now Eqs. (p^-|25|) reduce to 
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\ h n 36' u) + 



El 
IT 



El. (El 

U \2U 



V R' 

n h 

2V R 



m 



(- 

\2U 



El 

IT 



V R' 

h n — 

2V R 



E1(E1-Y1 

27 + 3/52 



/5 



0' 1 V^' + —^j 



V^ = o, 

(U' 



(28) 



(R'l±01J\ : 



'U' R' 

h n — 

2U R 



7 + 2/?- 
/3 



0' u 



(29) 



Using Eq. (|30|) , one can easily see that x can be evaluated from ip and -0, and that Eqs. (p8| , ^ 
become two second order coupled equations for Lp and ip only as follows 



r (r^ - k) if" + [{d + ly 



2 ~ 



2d ( r"' - — A;sinhV 



(J - 3) (r'^ + fcsinhV 



^2^d+2-4d/A 



A;) (r'' + fcsinh^/i 



1+4/A 



+ — (2rf2 + ((i + 3)(a2- 



2rf2 



sinh^ /i (^r*^ — k^ 



D-2 ^ 

Here the functions [/, and i? are substituted explicitly, and 

2 / , 2^2 



>V' + rfA;iy(^ = 0. (32) 



= J (r^ - A;) (r'^' + k sinh^ 



A 



sinhV(r'^-fc) +(iA:coshVr'^- (33) 



The question about whether the black p-branes are stable or unstable under linearized 
perturbations now becomes whether or not there exists some Kaluza-Klein mass parameter 
m for which the above coupled equations allow any regular solution outside the event horizon. 

Before going further, several points should be mentioned here. One may wonder if 
metric fluctuations alone with frozen dilaton perturbation could produce instability. When 
/ = 0, Eq. ( P^D gives —nx = fp + ip, and so from Eq. ( |5(]| ) one has an additional first 
order equation for if and ip. This raises the question about whether or not such equation is 
consistent with Eqs. (|3ll , |32D . It turns out that the inconsistency is proportional to sinh /i. 
We find therefore that metric fluctuations alone without dilaton perturbation cannot produce 
linearized instability when black branes are charged, i.e., /i 7^ 0. It is possible only for 
uncharged black branes, i.e., = 0. This property can also be expected by applying the 
same argument as in the work of Gregory and Lafiamme |jTO|. It is shown in Ref. [|10| that 
for small charge cases / = /o + ^{smh^ /i) where /o is independent of /i. Now it can be easily 
shown in Eq. (p2D that /o = 0. 

Secondly, Eqs. (^,^) are invariant under rescalings of m — > am, r„ — > ThIo., fi fi, 
and r ^ r/a for an arbitrary constant a. This scaling symmetry can be easily seen by 



defining dimensionless variables such as p = r/r^ and m 



mr. 



or by observing that 
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such rescalings with t ^ t/a and 2;' — z^/a simply result in an overall constant rescale 
of the metric under which field equations are not affected. This symmetry implies that 
the threshold mass must be inversely proportional to the horizon radius, i.e., m* ~ 1/t//. 
Accordingly, it suffices to study instability modes just for a single value of r„. 

Since perturbation equations above are coupled second order and linear in ip and 
there are four linearly independent mode solutions in general. The asymptotic solutions of 
Eqs. PT|P^) at the spatial infinity {i.e., r ~ 00) are given by 



d'+i (d—l)(d+l) d+3 

/■ 2 ip -i — 2 -|- ■ ■ 

8m 

d+3 (d + l)(d+5) d+3 

r 2 zp -f 2 + • • • 

8m 



(34) 
(35) 



up to overall arbitrary constants. By finding out asymptotic solutions in the vicinity of 
the event horizon as well, i.e., p = 1 + e, one can have four sets of mode solutions whose 
asymptotic behaviors are given by 

^ \ Aje-'^Pv.ip) + Bje^Pv+{p), ^" ^ \ A„e-'^^w_(p) + Bne^Pv+{p), ^ ' 

^ 1 Anre-'^'v^p) + B,,,e^Pv+{p), ~ \ A,,e-^Pv_{p) + B.^e'^Pv+ip), ^ ' 

and correspondingly 

/ e + ^{e^) / 1 + ^(cosh/i)8/^5 + ^(£2) 

"^'^ \Aje-'^Pu4p) + B,e^Pu+{p), ^ \ A„e-™^M_(p) + 5„e™''M+(p), ^ ' 

j elne-e + dieHns) f _ (coshpf Ine + ^^(e^) 

"^"^ ^ [ A,,,e-^Pu^{p) + B,,,e^Pu+{p), '^^^ ~ \ A,,e-^Pu.{p) + B.^-e'^Pu+ip), ^ ' 

where upper lines are for p ^ 1 and bottom lines for p ^ 00. All solutions are simply linear 
superpositions of these modes. Note that, for the mode IV , ipiv 7^ il^iv as one approaches 
to the horizon whereas for other three modes ip = ip. 

Now let us turn to the question of boundary conditions for solutions ip and if. Since 
they are linearized perturbations, they should be "small" at any positions outside the event 
horizon. Thus, at the spatial infinity we require both ip and ip should be the decaying 
one, that is, exponentially decreasing. In the vicinity of the horizon, the metric itself is 
regular in Kruskal coordinates if, and only if, the perturbation is bounded and ip = ip at the 



horizon [jT5[. All modes except for the mode IV above satisfy these conditions. However, it is 
very important to notice that this is not enough for the regularity of linearized perturbations. 
Note that the second derivative of the mode /// is proportional to \ne. Accordingly some 
curvature quantity, for instance, the perturbation of the Ricci scalar curvature associated 
with such mode, becomes singular at the horizon. In other words, the mode solution /// 
produces curvature singularity even if the metric perturbation itself is regular at the horizon. 
Mode solutions / and // do not produce any curvature singularity at the horizon. By 
requiring regular curvature at the horizon in addition, therefore, the boundary conditions 
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we impose for regular perturbations are that both ip and (p are hnear combinations of modes 
/ and // only near the horizon and that they should decay at the spatial infinity. 

For various values of parameters in Eqs. (^^2|), we have searched for regular solutions 
satisfying boundary conditions described above by using Mathematica. In more detail, we 
start from a regular solution at the horizon in the following form 

^P = C^pj + E^jj, ip = Cipj + Eifjj. (40) 

At the spatial infinity, they will become 

^P ~ (CA, + EAjj) e-'^'v^p) + {CBj + e™V(p), (41) 

if ~ {CAj + EAj,) e'"'Pu^{p) + (CB, + EBu) e^'^'uM- (42) 

By solving Eqs. (^,^) numerically, one can check if there exists any combination of con- 
stants, C, E, and fh such that coefficients of the exponentially growing parts for both and 
(/9 at the spatial infinity vanish in Eqs. ( |lT|j^ ), i.e., 

CBj + EB,j = 0, CBj + EB,j = 0. (43) 

However, we have used a different, but equivalent method that turned out to be more stable 
and efficient in actual numerical search. For instance, one just needs to check only for 
varying fh as shall be shown below. Defining 

P(m, /i, a, d, D) = BjBjj - BjjBj, (44) 



the existence of non-trivial constants C and E satisfying Eq. (|43| ) is guaranteed only if 
P = 0. Thus, having given the initial data as in Eqs. (^,^) near the horizon, we solve the 
coupled equations Eqs. (0,^) numerically, and evaluate the quantity P{fh, n, a, d, D) by 
using numerical values for ipuj and ipui at sufficiently large p. For given other parameters 
of black p-branes, we vary the Kaluza-Klein mass fh only and search for the value fh* at 
which P = 0. This can be achieved by finding out a fh* around which the function P above 
changes its sign. If there exists such m*, this mode is indeed the threshold unstable mode.i 

Once such fh* = fh*{p,a,d, D) is obtained numerically, the threshold mass m* = 
fh* /tjj = fh* /k^/'^ can be evaluated by fixing or k. Recall that a black p-brane in D 
dimensions is characterized by two physical quantities Q and M, or equivalently th and 

and one coupling parameter a. We would like to see how the stability behaves as the 
charge Q of a black brane increases when its mass density M is kept same.H We have chosen 
M such that threshold masses at /i = for the GL cases reproduce those in Gregory and 



Lafiamme's work [10,11|, i.e., tq = = 2. A non-extremality parameter can be defined as 



^In the actual numerical calculations, we used the sign change of ipn/ipi — ^n/y^i for p ^ 1, which 
becomes equivalent when Bjjj, Bjjj ^ 0. 

^In the Gregory and Lafiamme's work, they fix the size of the outer horizon tq in string frame 
and vary the inner horizon size r^. Here Tq = cosh^ p, and rf = rfj sinh^ p. 
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p = 4 

m* m* 




FIG. 1. Behavior of threshold masses for black 4-branes in D = 10 at various values of a with 
fixed mass density M = 2^ m* ~ LSSl/r^ ^ 0.791 with = 2 at /i = 0. ^ 0.818 (0.8184), 
0.881 (0.8814), and 1.125 (1.1254) for cases of a = 0, 1/2, and 1, respectively. Here critical values 
obtained from the GM conjecture are denoted in parentheses. On the righthand side the same 
data are plotted in terms of the non-extremality parameter q. n = 6 corresponds to g ~ 0.99996, 
0.99995, 0.99991 for a = 3/2, 2, 3. 



= _Q_ ^ J_Q_ ^ sinh2/i 
^~Q_„ v^M (J+ l)A/2(i + 2sinhV' 

Here Qmax is the maximum charge density allowed for a given mass density M. Since 
{d+ l)A/2d > 1, g is a monotonic function of fi and g = 1 is the extremal case. Notice that 
as or g increases, the horizon radius = k^^"^ decreases according to Eq. (|]) and becomes 
zero for the extremal brane when M is fixed. For black branes having different mass density 
M', one can get 

/ M \ ^/"^ 

m*{fi;M')=[^j^j m*{fx-M). (46) 

Behaviors of some of threshold masses we have obtained are illustrated in Fig. |I|. There 
we have considered black 4-branes in D = 10 and the mass density M = 2^. Our numerical 
results show that threshold masses at /i = are nonvanishing for any values of a, implying 
that all uncharged black 4-branes in Eq. (|^ are unstable under linearized perturbations 
associated with Kaluza-Klein massive modes with < m < m* as expected. In other 
words, instability occurs for small perturbations whose wavelengths in spatial worldvolume 
directions are in the range of A^,(= 27r/m*) < A < oo. The value of threshold masses at 
/i = can be seen to be independent of a as it should be so since the a-dependence of 
Eqs. (^,^) disappears as /i ^ 0. Numerically we find m* ~ 0.791 at /x = 0, which agrees 



with the numerical result obtained by Gregory and Laflamme [|10 . 

When black branes get charged, the background dilaton field becomes nontrivial {e.g., 
7^ 0) and could play some important role. As can be seen in Fig.|I|, the stability behaves 
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a| = 1/2 




FIG. 2. Behavior of threshold masses for black p-branes with fixed M in D = 10 in the theory 
of |a| = 1/2. At = 0, = 2 and m* = frf jr^ ~ 1.153, 1.044, 0.925, 0.791, 0.635, 0.440 for p = 1, 
• • •, 6. Critical values for transition points are =^ 0.413, 0.435, 0.497, 0.606, 0.773, > 0.9993 for 
p = l,---,6. 



very differently depending on the coupling parameter a between dilaton and gauge fields in 
Einstein frame Eq. (^. We summarize its behavior into three types, i.e., a < 3/2, a = 3/2 
and a > 3/2. 

In the work by Gregory and Lafiamme pil|JTT[ ], it has been shown that the black brane 
instability decreases monotonically as a black brane approaches to its extremal point. It 
is also noticed that numerical instability starts to occur near the extremal point, but the 
extremal brane was shown to be stable by analyzing its stability separately [|ll|]. The case 
for which numerical study was performed explicitly by Gregory and Lafiamme corresponds 
to the GL case in Fig. 0, i.e., a = —1/2 and p = 6. Our result confirms theirs up to ~ 4 or 
q ~ 0.9993. When /i > 4, numerical instability starts to occur as already observed in Refs. 
T^,TT]. For the GL case in Fig. ^ {i.e., a = 3/2 and p = 4), the stability behavior is similar 



at least up to /i ~ 4 or g ~ 0.99799. Our data for jj. > 4: are not reliable due to numerical 
instability. 

For cases of a = 0, 1/2, and 1, we find that threshold masses become zero at certain 
finite values of — 0.818, 0.881, and 1.125, respectively. Thus the Gregory-Lafiamme 
instability disappears as /i — ficr for such cases and so corresponding black 4-branes become 
in fact stable, at least under spherically symmetric linearized perturbations, for /icr < < 6, 
presumably all the way down to the extremal point {i.e., fi ~ oo). The same data are 



By decreasing instability, we mean that the threshold mass is decreasing. Thus, the instability 
actually shrinks in parameter range of and m. Consequently, smaller values of imply that 
instability modes grow less rapidly in time. So, in this sense we may say a black brane having 
larger charge is less unstable than those having smaller charges in the case considered by Gregory 
and Lafiamme. 
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plotted in terms of the non-extremality parameter q on the righthand side. Corresponding 
critical values for transitions in the stability are gcr — 0.593, 0.606, 0.674 for a = 0, 1/2, 1, 
respectively. This critical value gcr seems to increase monotonically from a minimum one 
for a = to unit for the GL case. Therefore one sees that there indeed exist some charged 
black branes which are stable even far from the extremality. 

For cases of a = 2 and 3, threshold masses decrease as black 4-branes get charged as 
before, but interestingly they start to increase again at around /i ~ 2 {i.e., q ~ 0.9972) and 
IX ~ 1.5 {i.e., q ~ 0.9950), respectively. This turning point becomes smaller as a increases. 
Thus it shows that the instability does not always reduce down as a black brane gains more 
charge. Such stability behavior in the presence of charge never has been expected in the 
literature as far as we know. Moreover, the threshold mass seems to diverge as the black 
brane approaches to the extremal point as can be seen in Fig. ^ up to g ~ 0.99995, 0.99991 
for a = 2, 3. As shown below, however, extremal black branes are stable in these cases as 
well. 

Now let us check the stability of extremal black branes. These cases are expected to 
be stable since they correspond to BPS ground states in string theory. Our numerical 
analysis cannot be directly used to show its stability since numerical instability occurs as 
the extremality is being approached. Instead, we treat the extremal case separately. As 
one can show explicitly, equations for linear perturbations in the extremal case are obtained 
simply by putting /c = and fi = oo with ke^^ = 4c fixed in Eqs. (|3ll , |32D , 



r 

(/ + c)V^" + 



^" + - ^ ( 1 + — ) = 0. (47) 

1 



r 



A 4 



,2/ 



m'(r- + c)(l + -^)^/^ 



i--^.W-^^V-^- (48) 




The boundary conditions are simply the regularity of and if), that are different from those 
of nonextremal cases, and the asymptotic solutions near the event horizon are also different 
and cannot be obtained simply by taking the extremal limit of those in Eqs. (p6|-p9|). This 
difference comes because we cannot ignore the terms including k at the near horizon as long 
as k is non-vanishing. Now it is straightforward to see that the only regular solution of 
Eq. ( ^71) is = since the coefficient of is negative definite. Similarly, the only regular 
solution for ip is zero. Therefore all extremal black branes considered in this paper are stable 
at least under s-wave fiuctuations, although the threshold mass seems to diverge for a > 3/2 
as the extremality is being approached. 

For black p-branes with different p (7^ 4), the basic stability behavior is essentially the 
same as the above. Only the values of threshold masses and critical ^ change slightly. Fig.|^ 
shows our results for various p-branes in the theory of \a\ = 1/2. The mass density is chosen 
to be M = (8 — p)2^~^, giving = 2 at /i = 0. If results for other values of a are added, 
one can see similar patterns as in Fig. ^ for each p-brane. Black p-branes with p = 1, ■ ■ • , 5 
becomes stable if /i > ficr whereas the instability of the 6-brane persists all the way down to 
at least = 4 {i.e., q = 0.9993). 

Now it will be very interesting to see how well the classical stability behavior of black 
branes described above agrees with that predicted by the local thermodynamic behavior 
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p=l 




oTs 1 iTs 2 275 3 Pl 

FIG. 3. Critical values of the parameter fi for various black p-branes in D = 10 at which the 
threshold mass vanishes m* = 0. The solid lines are obtained from the Gubser-Mitra conjecture 
and the black dots from our numerical results for several values of a. 

through the GM conjecture. First of all, let us consider critical values beyond that black 
branes become stable classically from being unstable. 



p 


1 


2 


3 


4 5 6 




(Num.) 
(GM) 


0.41^ 
0.41^ 


^ 0.549 
^6 0.5493 


0.695 
0.6954 


0.881 1.178 >4 
0.8814 1.1791 oo 



The table above shows the results for the case of \a\ = 1/2 in Fig.^. One can see that they 
are in very good agreement. Other cases we have checked are marked with black dots in 
Fig.^. Solid lines are obtained by using Eq. ([T7| ) in the GM conjecture. All critical values 
of n obtained numerically in our classical perturbation analysis agreed well with those in 
the GM conjecture. Based on such good agreement, this diagram shows how the classical 
stability under small perturbations will behave in general. That is, for a black brane with 
given dimension p for the spatial worldvolume, this brane will be unstable {i.e., /icr = oo) 
under small perturbations if |a| > Ocr = (-D — 3 — p)/^{D — 2)/2. If \a\ < Ocr, the brane is 
still unstable for /i < ficr{a,p, D), but becomes stable for /i > As can be seen in Fig.|^ 
for D = 10, critical values increase monotonically as \a\ or p increases. Critical values 
for various black brane solutions of the type II supergravity are marked in Fig. ^ explicitly. 



IV. DISCUSSION 



To conclude, we have investigated the stability of magnetically charged black brane 
solutions for low energy string theory in Eq. (|l|) under small perturbations. It turns out that 
all uncharged black branes in our consideration are unstable under linearized perturbations. 
When the black brane gets charged, however, the stability behavior depends on how strongly 
the n-form field couples to the dilaton field. In more detail, our results seem to show that, 
when the coupling is weak enough in the theory in Eq. (^ {e.g., \a\ < Qct), black branes 
become stable as they get charged enough {e.g., fi > ficr), even before they reach to the 
extremal point. When the coupling is strong enough {e.g., \a\ > a^), however, black brane 
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solutions of this theory are always unstable. Moreover the instability starts to increase 
again as the charge is larger than a certain value and seems to diverge near the extremality. 
For example, Fl, Dl, D2, D4 black branes of the type II supergravity could be stable 
classically for large charge, whereas NS5, D5 black brane solutions are always unstable all 
the way down to the extremal point. The case of a = is the boundary between these 
two categories, which is the case Gregory and Laflamme studied, and so the instability 
monotonically reduces down to zero at the extremal point. All extremal black branes are 
shown to be stable separately. It has also been shown that our results for the classical 
stability agree well with the qualitative behavior predicted by the local thermodynamic 
stability through the Gubser-Mitra conjecture. That is, the critical values for and ficr 
obtained numerically in our classical perturbation analysis agree very well with those in 
Eqs. ( p!6|JT7| ) for the sign change of the specific heat for the black brane being regarded as a 
thermal body. 

In the classical perturbation analysis above, we have considered only spherically symmet- 
ric perturbations. Thus, even if it turned out that there exists no instability mode for a black 
brane solution with fi > ficr in the analysis above, this does not necessarily mean such black 
brane is stable under small perturbations. There might exist some instability mode when we 
considered all non-s-wave perturbations as well. However, we give some evidence that the 
only possible instability mode comes from s-wave fluctuations as follows: Although it was in 



a different context, one can find in Ref. p5[ that the higher angular momentum fluctuations 



for Schwarzschild black brane backgrounds do not produce unstable modes. As /i — 0, since 
the black brane solutions considered in the present paper become Schwarzschild black branes 
and the metric perturbation equations become completely decoupled from others, one can 
easily see that the s-mode instability is the only instability for uncharged black branes in 
this paper. Now one can apply the same argument as in Ref. |TI|. When charge is added, 
we observed that there is a stabilizing influence for the s-wave perturbations if \a\ < a^r- 
Therefore, it is not expected that higher angular momentum modes exhibit instability in 
charged black branes since they do not give it even for the uncharged case. 

The spatial worldvolume is assumed to be non-compact in the description above. Con- 
sequently, the KK mass spectrum is continuous since there are translational symmetries in 
spatial worldvolume directions. What will happen on the stability behavior if the worldvol- 
ume is compactified with a scale LI In this case, the KK mass becomes discrete and has 
a minimum mass, e.g., mmin = 2'k/L. Thus, KK masses smaller than the minimum {i.e., 
< m < mmin) are not allowed due to the compactification. Accordingly, even if thresh- 
old masses obtained by solving perturbation equations for black branes are nonvanishing, 
such black branes are actually stable provided that m* < mmin, or equivalently that the 
compactification scale is small enough {i.e., L < 2n/m*) [P, p!0| , p!^] . 

Although critical values for transition points are in good agreement between the classical 
stability analysis and the local thermodynamic stability through the GM conjecture, there 
are some other aspects that might be in disagreement. For instance, one can see in cases 
of |a| < a^r that the magnitude of the specific heat increases whereas the black brane 
instability reduces as the charge increases with fixed mass density. Moreover, the specific 
heat is divergent at the transition point {i.e., n = Hcr)- In the classical perturbation analysis, 
however, there does not seem to exist any singular behavior around fi = /icr. For given 
/i in Fig.|I|, we found the threshold mass increases as a increases. One might wonder if 
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the magnitude of the specific heat also increases. One can see that it actually decreases 
monotonically, even for a fixed Q/M instead of a fixed ji. Maybe these discrepancies for some 
details in stability behavior simply indicate that one needs to find out some appropriately 
modified thermodynamic quantity, other than the specific heat, in order to compare other 
details between classical and thermodynamic stability analyses for black branes. 

Finally, it would be very interesting to understand why the instability starts to increase 
again as the extremality is being approached for black branes in theories of a > acr- Further 
work is required. 
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